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We discuss the static axially symmetric regular solutions, obtained recently in Einstein- Yang- Mills 
and Einstein- Yang-Mills-dilaton theory [1] . These asymptotically flat solutions are characterized by 
the winding number n > 1 and the node number k of the purely magnetic gauge field. The well- 
known spherically symmetric solutions have winding number n = 1. The axially symmetric solutions 
satisfy the same relations between the metric and the dilaton field as their spherically symmetric 
counterparts. Exhibiting a strong peak along the p-axis, the energy density of the matter fields of the 
axially symmetric solutions has a torus-like shape. For fixed winding number n with increasing node 
number k the solutions form sequences. The sequences of magnetically neutral non-abelian axially 
symmetric regular solutions with winding number n tend to magnetically charged abelian spherically 
symmetric limiting solutions, corresponding to "extremal" Einstein-Maxwell-dilaton solutions for 
finite values of 7 and to extremal Reissner-Nordstr0m solutions for 7 = 0, with n units of magnetic 
charge. 
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I. INTRODUCTION 



SU(2) Einstein- Yang-Mills (EYM) and Einstein- Yang-Mills-dilaton (EYMD) theory possess a rich variety of so- 
lutions. Recently, we have constructed static axially symmetric regular and black hole solutions in these theories 
1^,1). These solutions represent generalizations of the spherically symmetric regular and black hole solutions, known 
since several years All these solutions are characterized by two integers, the node number k of the gauge field 

function(s) and their winding number n with respect to the azimuthal angle <fi, since the fields wind n times around 
while 4> covers the full trigonometric circle once. The spherically symmetric solutions have winding number n = 1, a 
winding number n > 1 leads to axially symmetric solutions. 

All these EYM and EYMD solutions are asymptotically flat and have non-trivial magnetic gauge field configurations, 
but no global charge. To every regular spherically symmetric solution there exists a corresponding family of black hole 
solutions with regular event horizon ih > . The same holds true for the axially symmetric solutions J2| . Outside 
their event horizon these black hole solutions possess non-abelian hair. Therefore they represent counterexamples to 
the "no-hair" conjecture. The axially symmetric black hole solutions additionally contradict the expectation, that 
static black hole solutions are spherically symmetric. The spherically symmetric EYM and EYMD solutions are 
unstable and there is all reason to believe, that the axially symmetric solutions are unstable, too. 

For fixed n and increasing k the solutions form sequences, tending to limiting solutions. Although all solutions of a 
given sequence are neutral, the limiting solution is charged, and the winding number n of the solutions represents the 
charge of the limiting solution. The limiting solutions are spherically symmetric and abclian, representing Einstein- 
Maxwell-dilaton (EMD) [0 or Reissner-Nordstr0m (RN) solutions for EYMD and EYM theory, respectively. In 
particular, the limiting solutions of the sequences of regular solutions correspond to extremal EMD and RN solutions. 

The regular axially symmetric solutions have a torus-like shape. With the z-axis as symmetry axis, the energy den- 
sity has a strong peak along the p-axis and decreases monotonically along the z-axis jl]] . It has the same shape known 
from axially symmetric multimonopoles [RB, from axially symmetric multiskyrmions |]l0| and from multisphalerons 
|TT| . p^[ . The energy density of the axially symmetric black hole solutions looks torus-like only for small horizon radii, 
for larger horizon radii it changes shape [El . 

In EYMD theory the dilaton coupling constant 7 represents a parameter. In the limit 7^0 the dilaton decouples 
and EYM theory is obtained, for 7 = 1 contact with the low energy effective action of string theory is made, and in 
the limit 7—5-00 gravity decouples and Yang-Mills-dilaton (YMD) theory is obtained. The YMD limit is of particular 
interest, since YMD theory also possesses static spherically symmetric and axially symmetric solutions |l3|,|l2|]. In 
fact, the axially symmetric YMD solutions serve as starting configurations in our numerical construction of the regular 
EYMD and EYM solutions. 

This paper is the first of two papers on static axially symmetric solutions in EYM and EYMD theory and presents 
a detailed account of the regular solutions [Q . The second paper will give a detailed account of the axially symmetric 
black hole solutions [§,[l4| . In section II of this paper we present the action, the axially symmetric ansatz in isotropic 
spherical coordinates and the boundary conditions. Further we derive relations between the metric and the dilaton 
field and expressions for the mass. In section III we recall the spherically symmetric solutions, presenting them in 
isotropic coordinates. In section IV we discuss the properties of the axially symmetric regular EYM and EYMD 
solutions and the convergence of the sequences of regular non-abelian solutions to limiting abelian solutions. We 
present our conclusions in section V. In Appendix A we formulate the ansatz in cylindrical coordinates. We present 
the equations of motion in Appendix B and the expansion of the functions at the origin and at infinity in Appendix 
C. Finally, in Appendix D we discuss some technical aspects concerning the construction of the axially symmetric 
solutions. 



II. EINSTEIN- YANG-MILLS-DILATON EQUATIONS OF MOTION 



A. Einstein- Yang-Mills-dilaton action 

We consider the SU(2) Einstein- Yang-Mills-dilaton action 

Lm] 4~gd^x (1) 



with the matter Lagrangian 



(— 

\l61rG 



Lm = -\d^d^ - e 2K *±Tr( V") , (2) 
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the field strength tensor 

Fp, = d„A v - d v A ll + ie [Ap,A v ] , (3) 

the gauge field 

A, = \r a Al , (4) 

the dilaton field <&, and the Yang-Mills and dilaton coupling constants e and k, respectively. 
Variation of the action ([!]) with respect to the metric g^ v leads to the Einstein equations 

G>„ = R^ v — —g^i/R = SttGT^ (5) 

with stress-energy tensor 

.dL M 



Q gfJ .V 



= - \g^d a $d a <5> + 2e 2K *Tr(^ Q F^ ff ^ - ^g^F^F**) . (6) 

Variation with respect to the gauge field A^ and the dilaton field $ leads to the matter field equations. 



B. Static axially symmetric ansatz 

We here formulate the ansatz in terms of spherical coordinates Q . The equivalent formulation in terms of cylindrical 
coordinates jjj is given in Appendix A. 

Instead of the Schwarzschild-like coordinates, used for the spherically symmetric EYM and EYMD solutions |^-||] 
(see section III), we adopt isotropic coordinates to obtain the static axially symmetric solutions. In terms of the 
spherical coordinates r, 9 and (j> the isotropic metric reads ]l5| 

,9 „ , 9 171 , 9 mr 2 ,„9 If 2 sin 2 9 , , 

ds 2 = -fdt 2 + jdr 2 + —dO 2 + j # 2 , (7) 

where the metric functions /, m and I are only functions of the coordinates r and 9. The z-axis (9 = 0) represents 
the symmetry axis. Regularity on this axis requires [^6[ 

m\ e= o = l\e=o ■ (8) 
We take a purely magnetic gauge field, A = 0, and choose for the gauge field the ansatz (|Jl^,0,m,p| 



A^dx^ = ^- [r£ (Htdr + (1 - H 2 ) rd9) - n (t™# 3 + T%{l-H i ))r sin 0d0] . (9) 

Here the symbols t™, and denote the dot products of the cartesian vector of Pauli matrices, f — (t x ,t v ,t z ), 
with the spatial unit vectors 

e™ = (sin#cosn0, sin # sin n0, cos 9) , 
eg = (cos 6 cos n<j), cos sin n<j), — sin 0) , 

e? = (— sin 7i0, cosn0, 0) , (10) 

respectively. Since the fields wind n times around, while the azimuthal angle (f> covers the full trigonometric circle 
once, we refer to the integer n as the winding number of the solutions. The four gauge field functions Hi and the 
dilaton function $ depend only on the coordinates r and 9. For n = 1 and H\ = H3 = 0, H2 = H4 = w(r) and 
<I> = $(r), the spherically symmetric ansatz of ref. j^j is recovered. 

The ansatz (^) is axially symmetric in the sense, that a rotation around the z-axis can be compensated by a 
gauge rotation. The most general axially symmetric ansatz would contain nine gauge field functions (when Aq = 0). 
By requiring additional discrete symmetries this number is reduced. The ansatz (j^ respects the discrete mirror 
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symmetry M xz ® C .where the first factor represents reflection through the ir;2-plane and the second factor denotes 
charge conjugation |17| . |ll| , |l8| . 

The ansatz is form-invariant under the abelian gauge transformation J17|,pl],[l2[ 



(11) 



E/ = exp( -r£r(r,i 



The functions H\ and H 2 transform inhomogeneously under this gauge transformation, 

H X ^H X - rd r Y , 
H 2 ^H 2 + d e r , 



(12) 



like a 2-dimensional gauge field. The functions H 3 and H4 combine to form a scalar doubl et, (H 3 + ctgO^H^). The 
choice of an appropriate gauge leading to regular solutions presents a non-trivial problem |T^,|l9[. We therefore fix 
the gauge by choosing the same gauge condition as previously 1 17 T^JlJ,|l[ . In terms of the functions Hi it reads 



rd r H x - d e H 2 = 



(13) 



With the ansatz (@)-@ and the gauge condition (13) we obtain the set of EYMD field equations. These are given 
in Appendix B. Here we only present the energy density of the matter fields e = — Tq = — Lm 



1 o 



J_ 

2m 

n 2 

+ T 



2e 2 r 4 m [ m 



1 (rdrHi + deHi) 2 



[rd r H 3 - ffi£f 4 ) a + (rd r H 4 + H x (H 3 + ctgd)) 2 



+ (d H 3 -l + ctgeH 3 + H 2 H 4 f + (dgH 4 + ctge(H A -H 2 )-H 2 H 3 ) 2 ^ , (14) 
where the first gauge field term derives from F r g, the second and third derive from i 7 ^ and the fourth and fifth from 



C. Boundary conditions 



To obtain asymptotically flat solutions which are globally regular and possess the proper symmetries, we must 
impose certain boundary conditions [0,0] . Here we are looking for solutions with parity reflection symmetry. Therefore 
we need to consider the solutions only in the region < 9 < tt/2, imposing boundary conditions along the p- and 
^-axis (i.e. for — ir/2 and 9 — 0). The complete solutions then follow from the symmetries. In the following we 
discuss the boundary conditions for the regular solutions, which possess axial and parity reflection symmetry. The 
expansions of the functions at the origin (in powers of r) and at infinity (in powers of 1/r) are given in Appendix C. 

Boundary conditions at infinity 

Asymptotic flatness imposes for the metric functions of the regular solutions at infinity (r = 00) the boundary 
conditions 

f\r— 00 — 77l|r— 00 — ^\r— oo 1 • (1^) 

For the dilaton function we require 

$| r=oo = 0, (16) 

since any finite value of the dilaton field at infinity can always be transformed to zero via $ — > 3>(oo) , r — > re~ K *^ 00 - 1 . 
For magnetically neutral solutions, the gauge field functions Hi must satisfy 

H 2 \ r=00 = H4\ r=oc = ±1, -Hl| r =oo = #3|r=oo = . (17) 

Boundary conditions at the origin 

Requiring the solutions to be regular at the origin (r = 0) leads to the boundary conditions for the metric functions 

d r f\r=o = d r m\ r=0 = d r l\ r=0 = . (18) 
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Analogously, the dilaton function must satisfy 

d r $\ r=0 = , (19) 

and the set of gauge field functions Hi satisfies 

#2|r=0 = #4|r=0 = 1; #l|r=0 = #3^=0 = . (20) 

As for the spherically symmetric solutions , this choice of boundary conditions for the gauge field functions is 
always possible, because of the symmetry with respect to Hi — ► —Hi. Solutions with an even number of nodes then 
have 7?2(oo) = Hi(oo) = 1, whereas solutions with an odd number of nodes have ^2(00) = H^oo) = —1. 
Boundary conditions along the axes 

The boundary conditions along the p- and z-axis (9 — tt/2 and 8 = 0) are determined by the symmetries. The 
metric functions satisfy along the axes 



d e f\e=a = d s m\e=o = del\e=o = , 
<9e/|e=7r/2 = dem\ e=7T / 2 = d e l\e^/2 = , 

in addition to the condition (||) along the z-axis. Likewise the dilaton function satisfies 

<9e$|e=o = , 
d9$\e=n/2 = 

along the axes. For the gauge field functions Hi symmetry considerations lead to the boundary conditions 
Hi\g =0 = H 3 \e =0 = , 9ei?2|e=o = deH^g^o = , 

#l|e=7r/2 = #3|e=7r/2 = > deH 2 \e=Tr/2 = deHi\g =1T /2 = 

along the axes. 



(21) 



(22) 



(23) 



D. Dimensionless Quantities 

It is convenient to introduce dimensionless quantities. We introduce the dimensionless coordinate x, 



x = -^==r , (24) 



the dimensionless dilaton function ip, 



ip = Vi^G$ , (25) 
and the dimensionless dilaton coupling constant 7, 

7 =vir- (26) 

The value 7 = 1 corresponds to string theory. 

The dimensionless mass /1 is related to the mass M via 

" = ^ M ' (27) 
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E. Relations between metric and dilaton field 



We now derive two relations between metric and dilaton field for static axially symmetric solutions of EYMD theory 
with purely magnetic gauge fields and a general gauge group. These relations are known to hold for the spherically 
symmetric solutions pO^ 

Let us return to eq (|) and introduce the gauge field Lagrangian La 



La 



-Tr(F^F^ 



From the action (|l|) we obtain the equation of motion for the dilaton field 



(28) 



(29) 



Our aim now is to replace the matter field functions on the r.h.s. of this equation by metric functions alone. Therefore 
we consider the matter lagrangian 



L 



1 



M 



We replace Lm on the l.h.s. with help of the (00) component of the Einstein equations 

Goo = 87rGT o = S-rrGgooLM , 
and we replace the dilaton term on the r.h.s. with help of the contracted Einstein equations, 

R = 



inGTj* 



InGd^d^® 



(30) 



(For purely magnetic gauge fields the contracted gauge field stress energy tensor vanishes.) This gives for the matter 
part on the r.h.s. of eq. (E9h 



2e^L A = 2 ff uu T 00 - g^T^ 



1 



8ttG 



(2.9 00 G 



oo 



R) 



l 



4ttG 



9 Moo 



and the dilaton field equation becomes 

d u 



-gd^) 



AttG 



-99°°Roo 



By evaluating i?oo explicitly for the metric (0), we find 



-5.g°°i?oo = ^(J—gdnnf) 



The dilaton equation then becomes 



0, 



8ttG 



In/ 



= 



The derivation of the relations between metric and dilaton field is based on this equation. 
Abbreviating the expression in square brackets in eq. ( [34]) by B, 



B = $ 



inG 



In/ , 



(31) 



(32) 



(33) 



(34) 



(35) 



we now show that the regular static axially symmetric solutions satisfy B = 0. Therefore we multiply eq. ( |34"| ) by B 
and integrate over r and 9. When integrating by parts, the surface terms vanish, (the boundary conditions eq. ( |15| ) 
and ( |l^ ) give B(oo) = 0,) and the integral becomes 



/■OO /*7T 

/ / (r 2 {d r Bf + {d e Bf) sin 9Vldrd9 = 
Jo Jo 



(36) 



Since the integrand is positive this implies that B is constant. With B(oo) = this implies B = 0, yielding the first 
relation between metric and dilaton field, 



G 



Using dimcnsionlcss quantities, this relation reads 



8ttG 



In/ 



7 7 

representing the same relation obtained in the case of regular spherically symmetric solutions [ p0[ 
The second relation is either obtained by integrating eq. ([54|) over r and 8, or directly from eq. 

lim r 2 d r <& = — ^— lim r 2 <9 r / , 

r — >oo $7rG r — >°° 



since /(oo) = 1. Changing to dimensionless quantities we find 

1 ] 

2 x 



lim x 2 d x tp = — lim x 2 d x f . 



The l.h.s. of eq. ( f40| ) defines the dilaton charge 



D = lim x 2 d x ip 



The r.h.s. of eq. (40) is proportional to the dimensionless mass fi 

\i = - lim x 2 d x f , 

Z X — >OG 

as shown below. Together this yields the second relation 

which also holds for regular spherically symmetric solutions pQj. 



(37) 
(38) 

|), which yields 
(39) 

(40) 
(41) 

(42) 
(43) 



F. Mass 



The mass M of the static axially symmetric solutions can be obtained directly from the total energy-momentum 
"tensor" of matter and gravitation pif 

M = J T ao d 3 r . (44) 

Adopting the procedure of we start from the metric g^ in isotropic coordinates, eq. (Q). We then change 
to quasi-Minkowskian coordinates, in the sense that the metric g^ approaches the Minkowski metric rj^y at great 
distances, i.e. we introduce h^ LV 

9nu = + h^u (45) 
such that hfj, v vanishes at infinity. Identifying from the metric equation, eq. (^), and employing the coordinates 



x 1 = r sin# cos , x 2 — r sin 9 sin , x 3 =rcos9, (46) 



we obtain 



\ 



V 



Of — j^j sin cos 




fcos 2 (j)+ jsin 2 0- 1 - i ) sin0cos0 






f sin 2 + j cos 2 0-1 

f-l 



(47) 



Indeed, h^ v vanishes at infinity. 

The exact Einstein equations can be written as 
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11 fiv 9 '/(U!/- f T- A 



a - ~8^G (T M „ + V) - (48) 

where R is the part of the Ricci tensor linear in /i^. Interpreting the quantity 

r af3 = rf*»tf v (T M „ + t M „) (49) 

as the total energy-momentum "tensor" of matter and gravitation, the mass is given as the volume integral over r 00 . 
Transforming this integral to a surface integral yields for the mass J2lj 



M —i6^J{^-^t)i rsMd ^ ^ 

where x 1 are the coordinates, eq. (^), and the integral is taken over a large sphere of radius r — > oo. Evaluating the 
integrand, we obtain 

1 f f d f m + l\ m — l\ 9 
= - 16^ / [d-r {—) + "TtJ ^ Sin ^ • < 51 > 

Since the term in square brackets is independent of the angle for r — > oo, the integration is trivial, yielding 

1 f d f m + l\ m — l\ n 1 o„ , „ 

where the last step makes use of the asymptotic behaviour of the metric functions (see Appendix C) 

/ = l + - + o(-V) , m = l + o(V) > ' = 1 + f^ • (53) 

In dimensionless quantities, eq. ( f52] ) yields the desired relation for the mass, eq. (|42|). 

We next show, that the mass of the regular axially symmetric solutions can also be obtained from 

M = — / / (2T ° - T/) drddd0 . (54) 

Jo Jo Jo 

We start from the dilaton field equation ( |32| ) 

9, (V=S***) = ^a^^f (55) 

and replace the r.h.s with help of eq. (^Tj) 

8, (V^gd^) = -K^—g (2T ° - T/) . (56) 
Integrating both sides over r, 9 and 0, we obtain 

4tt lim r 2 <9 r $ = -k / / / (2T ° - T ' ») drd8d(f> = kM , (57) 

r ^°° Jo Jo Jo 



where the last equality makes use of relation (|39|) together with eq. ( |52[ ) . 

III. SPHERICALLY SYMMETRIC SOLUTIONS 

The spherically symmetric EYM and EYMD solutions were obtained previously in Schwarzschild-like coordinates 
with metric ||f ,§§ 

ds 2 = -A 2 Ndt 2 + j^dr 2 + f 2 (d6 2 + sin 2 8d4> 2 ) , (58) 
where the metric functions A and N are only functions of the radial coordinate f, and 

m = l _ ^ . (59) 
r 

Since we here construct the static axially symmetric solutions in isotropic coordinates, we briefly review the spherical 
solutions in isotropic coordinates, derive the coordinate transformation between both sets of coordinates and discuss 
the limiting solutions. 
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A. Coordinate transformation 



By requiring I — m and the metric functions / and m to be only functions of the coordinate r, the axially symmetric 
isotropic metric (|t]) reduces to the spherically symmetric isotropic metric 



Comparison with the metric 



and, with (plj), 



ds 2 = -fdt 2 
3) yields 



771 (A 2 

j(dr 



sm 



6d</> 2 )) 



m(r) 

7M 



dr 
r 



dr 



(60) 



(61) 



(62) 



The function N(f) (or equivalently the mass function m(f )) is only known numerically. Therefore the coordinate 
function r(r) can only be obtained numerically from eq. (p2|). To obtain the coordinate transformation, we switch to 
dimcnsionless coordinates x and x and introduce the function (3, 



Rewriting eq. (|62[) in terms of (3 



we find 



d\n(3 
dx 



x 
x 



N 



j3(x) = exp 



1 dx' 



with 



N{x) 



1 



2fi(x) 



(63) 



(64) 



(65) 



(66) 



The integrand is well behaved at the origin, since n(x) ~ i 3 , and the integration constant is adjusted such that at 
infinity [3=1, i.e. x = x. 

Fig. 1 demonstrates the coordinate transformation for the spherically symmetric EYMD solutions, for 7=1 and 
k = 1 — 4. Figs. 2a, b show the metric functions / and m, and Figs. 3a, b show the gauge field function w and the 
dilaton function ip. These solutions satisfy relation (p^), ip = (7/2)ln/ and relation (p3|), D = 7/1, where /i = /x(oo) 
represents the dimensionless mass. 

Fig. 4 demonstrates the coordinate transformation for the spherically symmetric EYM solutions, for k = 1 — 4. 
Figs. 5a, b show the metric functions / and m. 



B. Limiting solutions 

For fixed dilaton coupling constant 7, the sequence of neutral regular spherically symmetric EYMD solutions 
converges to a limiting solution J^,^]. This limiting solution is the "extremal" EMD solution Q| with horizon x = 0, 
charge P = 1, mass 



/i = P/^l+7 2 (67) 
and dilaton charge _D = 7/i (|38|). For 7 = 1 the coordinate transformation for the limiting solution 



-\= (V-P 2 + 2z 2 - P) (68) 
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is also shown in Fig. 1, the metric functions / and m of the limiting solution 



d S 2 = -[ —± ) dt 2 + f ^ P + X ) ^ + ^ (^2 + ^2 ^2)) (6Q) 



are also shown in Figs. 2, and the dilaton function if of the limiting solution, 

e 2 ^ = , (70) 

V2P + x 

is also shown in Fig. 3. The gauge field function of the limiting solution is trivial, = 0. The limiting functions for 
/ and ip are reached uniformly. In contrast, the function m approaches the limiting function m,oo = 1 non-uniformly, 
since it tends to a constant m(0) ^ 1 at x = 0. Likewise the function w approaches the limiting function Woo = 
non-uniformly, since the boundary conditions require w(0) ^ as well as w(oo) ^ 0. 

The sequence of spherically symmetric regular EYM solutions also approaches a limiting solution [Q. In 
Schwarzschild-like coordinates this limiting solution consists of two parts. A non-trivial part for x < 1, representing 
an oscillating solution, and a simple part for x > 1, representing the exterior of an extremal RN solution with mass 
[1=1, horizon radius xh = 1 and charge P = 1. 

From Fig. 4 we observe, that with increasing k the interval < x < 1 is mapped to an ever smaller interval 
< x < x m , with x m tending to zero for k — > oo. Therefore the limiting solution becomes rather simple in isotropic 
coordinates. The non-trivial part of the solution in Schwarzschild-like coordinates is mapped to zero, and the limiting 
solution corresponds to the exterior of an extremal RN solution with mass [i = 1, horizon ih = and charge P = 1. 
The coordinate transformation for this limiting RN solution, 

x=x-P, (71) 

is also shown in Fig. 4, and the metric functions /oo and moo, given by the metric 



ds 2 



(-p^-) 2 d e + (^1) 2 (dx 2 + x 2 (dO 2 + sin 2 Odcf)) (72) 



(i.e. moo = 1), are also shown in Figs. 5. Again, the gauge field function of the limiting solution is trivial, w c 
and the functions m and w approach the limiting functions non-uniformly. 



IV. AXIALLY SYMMETRIC SOLUTIONS 



Beginning with a short discription of the numerical technique, we here give a detailed description of the numerically 
constructed axially symmetric regular solutions. We discuss the properties of these solutions and the convergence of 
the various sequences of solutions to limiting solutions. Some technical aspects concerning the construction of the 
axially symmetric solutions are discussed in Appendix D. 



A. Numerical technique 

Subject to the above boundary conditions, we solve the system of coupled non-linear partial differential equations 
numerically. We employ the radial coordinate 

(73) 



1 + x 



instead of x, to map spatial infinity to the finite value x = 1 (see also Appendix B). The numerical calculations are 
performed with help of the program FIDISOL, which is extensively documented in p3j. The equations are discretized 
on a non-equidistant grid in x and 0. Typical grids used have sizes 150 x 30, covering the integration region < x < 1 
and < 9 < it/2. 

The numerical method is based on the Newton-Raphson method, an iterative procedure to find a good approximation 
to the exact solution. Let us put the partial differential equations into the form P(u) — 0, where u denotes the unknown 
functions (and their derivatives). For an approximate solution u^, P(u^) does not vanish. If we could find a small 
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correction Am, such that + Au is the exact solution, P(u^ + Au) — should hold. Approximately the expansion 
in Am gives 

BP 

= P{u {1) + Au) « P(u {iy ) + |-(« (1) )Am . 

The equation P(u^) — — ^(u^')Au can be solved to determine the correction Au^ = Au. = + Au^ 
will not be the exact solution but an improved approximate solution. Repeating the calculations iteratively, the 
approximate solutions will converge to the exact solutuion, provided the initial guess solution is close enough to the 
exact solution. The iteration stops after i steps if the Newton residual P(u^') is smaller than a prescribed tolerance. 
Therefore it is essential to have a good first guess, to start the iteration procedure. Our strategy therefore is to use 
a known solution as guess and then vary some parameter to produce the next solution. To construct regular axially 
symmetric EYMD and EYM solutions, we have used the known regular axially symmetric YMD solutions JlJ] as 
starting solutions, corresponding to the limit 7 — > 00. Beginning with large values of the parameter 7, 7 is slowly 
decreased, until EYMD solutions with 7 = 1 and EYM solutions (7 = 0) are obtained. A second procedure, also 
employed, is to start from the spherically symmetric regular solutions with n — 1, and then increase the 'parameter' 
n slowly, to obtain the desired EYMD and EYM solutions at integer values of n. 

For a numerical solution it is important to have information about its quality, i. e. to have an error estimate. The 
error originates from the discretization of the system of partial differential equations. It depends on the number 
of gridpoints and on the order of consistency of the differential formulae for the derivatives. FIDISOL provides an 
error estimate for each unknown function, corresponding to the maximum of the discretization error divided by the 
maximum of the function. For the solutions presented here the estimations of the relative error for the functions are 
on the order of 10 -3 and 10" 2 for k < 4 and k — 4, respectively. 

B. Properties of the solutions 

For the static axially symmetric solutions constructed here, the functions depend on two variables, the radial 
coordinate x and the angle 0. In the following we either present the functions in three-dimensional plots, employing 
the compactified coordinates p — xsinO and z — xcos0 (with x, given in eq. (|73|)), or we present the functions in 
two-dimensional plots exhibiting the x-dependence for three fixed angles, = 0, 6 — 7r/4 and = ir/2. While the 
three-dimensional plots give a good qualitative picture of the solutions, the two-dimensional plots are better for a 
quantitative analysis of the solutions and in particular for a comparison of solutions with differing values of n, k or 7. 

We begin with a description of the lowest axially symmetric regular solution of EYMD and EYM theory. This 
solution has winding number n = 2 and node number k = 1. The EYMD solution for 7 = 1 is shown in Figs. 6-9. 
Figs. 6a-c show the metric functions /, m and I. The functions m and I have a rather similar shape. Figs. 7a-d show 
the gauge field functions H1-H4. The single (non-trivial) node of the functions H2, H3 and H4 is seen in Fig. 7e. The 
lines representing the position of the node of H2 and H4 are close together. The non-trivial node of the function H3 
is located farther outwards. Interestingly, the nodal lines of the functions H2-H4 are almost spherical. The function 
Hi does not possess a non-trivial node. Fig. 8 shows the dilaton function ip, which is related to the metric function 
/ via (|38|). Fig. 9 presents the energy density of the matter fields e, showing a pronounced peak along the p-axis 
and decreasing monotonically along the z-axis. Thus equal density contours reveal a torus- like shape of the solutions. 
The maximum of the energy density e max and the location p max of the maximum of the energy density are shown in 
Table 1 along with the mass of the solution. 

The corresponding EYM solution looks very similar to the EYMD solution. We therefore exhibit only the energy 
density of the matter fields for this solution in Fig. 10. The mass of the solution is shown in Table 2 along with the 
maximum of the energy density e max and the location p max of the maximum of the energy density. 

Let us now consider the lowest solutions for higher winding number. To see the change of the functions for node 
number k = 1 with increasing winding number n, we exhibit in Figs. 11-14 the EYMD solutions for 7=1 with n = 1, 
2 and 4. Figs, lla-c show the metric functions, Figs. 12a-d the gauge field functions, Fig. 13 the dilaton function and 
Fig. 14 the energy density of the matter fields. Obviously, the angle-dependence of the metric and matter functions 
increases strongly with n. The location of the biggest angular splitting of the metric and the matter field functions 
moves further outward with increasing n along with the location of the nodes of the gauge field functions H2-H4. 
At the same time the peak of the energy density along the p-axis shifts outward with increasing n and decreases in 
height. At the origin the values of the metric functions increase with n, while the central energy density decreases. 
The nodal line of the functions H2-H4 remains almost spherical with increasing n, and Hi remains nodeless. (Note, 
that Hi and H3 are zero on the axes in Figs. 12a, c.) 



11 



By decreasing 7 from one to zero, we obtain the corresponding EYM solutions, shown in Figs. 15-17. Qualitatively 
the EYM solutions look like their EYMD counterparts. But there are a number of quantitative differences. The 
metric functions of the EYM solutions are considerably smaller at the origin, the gauge field functions have their 
peaks and nodes shifted inwards, and the energy density of the matter fields has slightly higher peaks, which are also 
shifted inwards, as compared to the EYMD solutions. The masses increase with decreasing 7, as seen in Tables 1-3. 

We now consider the radially excited solutions. To see the change of the functions for fixed winding number n with 
increasing node number k, we exhibit in Figs. 18-21 the EYMD solutions for 7=1 winding number n = 2 and node 
number k = 1 — 4. Figs. 18a-c show the metric functions, Figs. 19a-d the gauge field functions, Fig. 20 the dilaton 
function and Fig. 21 the energy density of the matter fields. The magnitude of the angular dependence of the metric 
and matter functions stays roughly the same with increasing k in these figures. The peak of the energy density of 
the matter fields moves inward with increasing k and at the same time increases strongly in height, along with the 
central density. The gauge field functions H2 and H4 have precisely k nodes, the gauge field function Hi has k — 1 
nodes, and the gauge field function H3 has one node. 

Again, the corresponding EYM solutions are similar to these EYMD solutions. We exhibit only the metric functions 
in Figs. 22 and the energy density of the matter fields in Fig. 23. Again, the metric functions of the EYM solutions 
are considerably smaller at the origin, and the gauge field functions have their peaks and nodes shifted inwards, as 
compared to the EYMD solutions. But the energy density of the matter fields increases only moderately with k, 
indicating a finite limit for the energy density of the EYM solutions. In contrast, the dramatic increase of the energy 
density of the matter fields of the EYMD solutions with increasing k indicates a diverging limit. 

The masses of the EYMD solutions with 7=1 are shown in Table 1, and the masses of the EYM solutions are 
shown in Table 2. The EYMD solutions for other values of 7 are similar to these solutions. For comparison, Table 3 
presents the masses of the solutions with n = 3 and k = 1 — 4 for several values of 7. In the limit 7 — > the EYM 
solutions are obtained, and in the limit 7 — > 00 the corresponding YMD solutions [p]Jl^] . 

C. Limiting solutions 

Let us now consider the limiting solutions of the sequences of axially symmetric solutions with fixed winding number 
n, dilaton coupling constant 7 and increasing node number k. For n > 1 the pattern of convergence is similar to the 
one observed in the spherically symmetric case, n = 1, considered in section III. 

For fixed winding number n and dilaton coupling constant 7, we consider the solutions with k = 1 — 4, shown in 
Figs. 18-23 and Tables 1-3, and extrapolate to the limit k — > 00, to obtain the limiting solution, also shown. For finite 
7 this limiting solution is an "extremal" EMD solution with mass fioo = n/ 'yl + 7 2 (|6^), charge P = n and metric 
(|69|). For 7 = it represents the exterior solution of an extremal RN solution with mass fi^ — n, horizon xjj = 0, 
charge P — n and metric (|7l|). So we find, that this limiting solution is spherically symmetric, although the solutions 
of the sequence are axially symmetric, and that this limiting solution carries charge n, although the solutions of the 
sequence are magnetically neutral. 

The limiting values for the mass represent upper bounds for the sequences. The convergence of the mass fik to the 
mass of the corresponding limiting solution [loo is exponential. This is demonstrated in Figs. 24, where we show the 
deviation of the mass of the fc-th solution from the mass of the limiting solution, A^., 

A/i fc = ^00(00) - ^^(oo) , (74) 

for fixed winding number n, n = 1 — 4, and 7=1 (Fig. 24a) as well as 7 = (Fig. 24b). As required for an exponential 
convergence, the logarithm of the deviation A^ik is a straight line as a function of fc. The slope of this straight line 
decreases with increasing n, thus the larger n the slower is the convergence to the limiting solution. 

The convergence of the functions of the axially symmetric solutions is less rapid than the convergence of the mass. 
The metric function /, the gauge field functions H\ and H3 and the dilaton function tp converge uniformly to the 
corresponding functions f OQl i?i i0 o, ^3,00 and t^oo of the limiting solution. For instance we observe from Fig. 20, that 
the dilaton function (p deviates from ipoa only in an inner region, which decreases exponentially with k. The metric 
functions m and I also deviate from the limiting functions = 1 and = 1 only in an exponentially decreasing 
inner region. But the metric functions m and I as well as the gauge field functions H2 and H4 converge non-uniformly 
to the corresponding limiting functions. The metric functions m and I tend to constants m(0) ^ 1 and 1(0) ^ 1 at 
the origin, while the gauge field functions H2 and H4 differ from the value zero of the limiting solution both at the 
origin and at infinity because of their boundary conditions. 
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V. CONCLUSIONS 



We have given strong numerical evidence, that both EYMD and EYM theory possess sequences of regular static 
axially symmetric solutions, in addition to the known static spherically symmetric solutions. These sequences are 
characterized by the winding number n > 1, and the solutions within each sequence by the node number k. (For 
n = 1 the spherically symmetric solutions are recovered.) To rigorously establish the existence of these solutions a 
formal existence proof is desirable, analogous to the proof given for the spherically symmetric solutions 24 1. 

For fixed n and 7 > (7 = 0) the sequences of regular axially symmetric solutions tend to limiting solutions. These 
are the "extremal" EMD solutions j?J (the exterior solutions of the extremal RN solutions) with magnetic charge n 
and the same 7. This generalizes the corresponding observation for the sequences of spherically symmetric solutions 
P)P0|. For the YMD solutions the limiting solutions are magnetic monopoles with n units of charge fl^| . 

The regular axially symmetric solutions have a torus- like shape. Otherwise, many properties of the axially symmetric 
solutions are similar to those of their spherically symmetric counterparts. In particular, there is all reason to believe, 
that these static regular axially symmetric EYMD and EYM solutions are unstable. Since we can also associate 
the Chern-Simons number -Yes = n /2 |]ll| , |l8| (for odd k |p5| ) with these solutions, we interpret them as gravitating 
multisphalerons. 

Having constructed the axially symmetric solutions in EYMD and EYM theory, it appears straightforward to 
construct analogous solutions in other non-abelian theories. For instance, we expect to find gravitating axially 
symmetric multimonopoles in SU(2) Einstein- Yang-Mills-Higgs (EYMH) theory with a Higgs triplet, gravitating 
axially symmetric multisphalerons in SU(2) EYMH theory with a Higgs doublet, or gravitating axially symmetric 
multiskyrmions in SU(2) Einstein-Skyrme (ES) theory. In particular we expect, that the lowest gravitating axially 
symmetric n = 2 multimonopole and multiskyrmion solutions are stable. 

We expect, that the axially symmetric solutions represent only the simplest type of non-spherical gravitating regular 
solutions in EYMD and EYM theory as well as in other non-abelian theories. We conjecture, that there are gravitating 
regular solutions with much more complex shapes and only discrete symmetries left. Such solutions arise in flat space 
for instance in SU(2) Skyrme theory |2(|^7j and in SU(2) Yang-Mills-Higgs theory pj| , representing multiskyrmions 
with higher baryon number and multimonopoles with higher charge, respectively. 

But EYMD and EYM theory also possess black hole solutions. The non-abelian spherically symmetric black hole 
solutions may be regarded as black holes inside sphalerons plfl . The axially symmetric black hole solutions in EYM 
and EYMD theory corresponding to black holes inside multisphalerons, will be discussed in detail in ]l4| . We 
expect analogous black hole solutions in EYMH and ES theory. Indeed, a perturbative calculation has revealed not 
only axially symmetric EYMH black holes, but also black holes without rotational symmetry p9[ . 

We would like to thank the RRZN in Hannover for computing time. 
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VI. APPENDIX A 



We here present the ansatz, the gauge condition and the boundary conditions in cylindrical coordinates, as employed 
in 111]. We also show the relation of the functions Hi, used in this paper, to the functions i 7 !;, employed in 0. 



A. Ansatz in cylindrical coordinates 



In terms of the cylindrical coordinates p and z the metric in isotropic coordinates reads 



ds 1 



-fdt 2 + - (dp 2 + dz 2 ) + jp 2 d(f> 2 



(75) 



where the metric functions /, m and I arc only functions of p and z. 
The gauge field is parametrized by |8| Ju)Jrl| |l| 



A*^" = - [t£ (w\dp + w\dz) + {r£w\ + T z wfj po 



(76) 



where t? , r™ and t z denote the dot products of the cartesian vector of Pauli matrices with the spatial unit vectors 



e'£ = (— sinn0, cosn</>, 0) , = (cos ncj), sin ruj), 0) 



(0,0,1) , 



(77) 



respectively. The four functions and the dilaton function <£> depend only on the coordinates p and z. 
With respect to the gauge transformation (111]), 



C/ = exp( -t$T(j>,z) 



(78) 



the functions (pw 3 , pw 2 —n/e) transform like a scalar doublet, and the functions (wf, w 2 ) transform like a 2-dimensional 
gauge field. The gauge condition reads |p7| , |ll| , ^2| Jl| 

d p w\ + d z wl = . (79) 

With the ansatz (|7^)-(76) and the gauge condition (|79| ) the set of EYMD field equations is obtained. The energy 
density of the matter fields e = — Tq = —Lm reads 



2k<& f 2 

2ml 



dpwl + - - — - ew'{w 2 

2 



d p wl 



d z w 3 



3 \ z 

i TIW2 3 2 \ / o 2 3 1\* 

3 H ew 2 w 3 J + {OzW 3 + ew 2 w 3 ) 



„3„„1 



P 3 



(80) 



Transforming to spherical coordinates r and 9 (p = rsin^, z = rcos9), the gauge field functions Fi{r, 
introduced via 



arc 



— (l-Fi)cos6, w\ 

er 

- — (1 - F 3 )cos6> , wl 



— (1-F 2 )sin0 , 
er 

TL 

— (1 - F 4 )sin6> . 



(81) 



The spherically symmetric ansatz of ref. || is recovered for Fi — F 2 — F 3 — F4 — w(r), <I> = $(r) and n = 1. 
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B. Boundary conditions 

At infinity the boundary conditions for the metric functions and the dilaton functions are given by eqs. (|l5|)-([i"6|), 
whereas the gauge field functions Fi satisfy 

•Pl|r=oo = ^2|j-=oo = ■P3|r=oo = ^*4|r=oo = ±1 ■ (82) 

At the origin the boundary conditions for the metric functions and the dilaton functions are given by eqs. (p^)-(^9|), 
and the set of gauge field functions Fi satisfies 

Fi\ r= Q = F-2\ r =0 = ^3^=0 = -^4^=0 = 1 ■ (83) 

Alon g th e p- and z-axis, the boundary conditions for the metric functions and the dilaton function are given by 
eqs. (p2|) , and for the gauge field functions Fi by 

dgFi\g=o = dgF 2 \e=o = dgF 3 \g= = dgF 4 \g= a = , 

d$F 1 \ e=7I /2 = deF 2 \ s= T I /2 = dgF 3 \ g=7I /2 = 0e-p4|e=7r/a = 
with the exception of F 2 (r, 9) for n = 2, which has 

dgF 2 {r,9)\ e=0 ^0 (tarn =2). (85) 



C. Relation between Hi and Fi 

The functions Fi are related to the functions Hi (r, 9) via 

Hi = sm0cos0(F 2 - Fi) , 
H 3 = sin 9 cos 6 (F 4 - F 3 ) , 
H 2 = cos 2 6F X + sin 2 9F 2 , 

H A = cos 2 8F 3 + sin 2 9F 4 . (86) 

Note, that the boundary conditions on the z-axis for the functions Fi, eqs . (p4|) , imply the boundary conditions for 
the functions Hi, eqs. (p3|) , but not vice versa. The boundary conditions fl23|) allow for dgF 2 (r,9)\g=o ^ i. e. they 
properly incorporate the case n = 2 with condition (|S5|). 

VII. APPENDIX B 

We here derive the equations of motion in isotropic spherical coordinates. 

A. Tensors F^, T MlM 

We expand the field strength tensor with respect to the Pauli matrices r™, (A = r, 9, 4>, see eq. (|lC|)) 

F - p W ZA. 



Inserting ansatz (|J) for the gauge field, we obtain the non-vanishing coefficients F$ 



F, 



F, 



(r 



F 



(0 



F 



(0 



--(H h g+rH 2 , r ) , 
r 

-n—irH^-HrH* 



Dill \J 

n (rH 4 r + HiH 3 + cot 6H X ) , 

r 

-n sin 9 (H 3 . g - 1 + H 2 H 4 + cot 8H 3 ) , 
n sin 9 (H± fi - H 2 H 3 - cot 9 (H 2 - H A )) , 



(87) 
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and Fjw' — —Fi>n - It is convenient to define 

F r$ = (^rV) (^V) ' 
F H = ( F w) + (^eT ) ' 



where the second term in the definition of F^ g represents the gauge fixing term. 
With the ansatz for the metric (Q) we obtain the Lagrange densities 

L * = -T-(**r + ?s**e) . 

2 m \ r z J 



and the non- vanishing components of the stress energy tensor, 



T 



T, 



1 



2to 



- $i - — $ 



2 | g2ft$ 



2 | g 2K* 



/ 7Tl2 

"2 t r8 

t; — iVs + 



/ 



1 



> • J nS F r<t> + ~^ F 8<b 

r A sin 0/ r z v 

1 



2 

r 2 sin 
2m 



2 | g 2K* 



/ F 2 



f C^r<£ ~ ~ F S4>) 



r2 sin 2 0/ ^ r0 r 2 



_<J> Z <T> 2 



r Tn 

,2k$ 



r 2 sin 2 61 



~2 b rd 



f 1 



r 2 sin 2 6i^ + r- F 



2k$ 



/ 

r 2 sin 2 0Z 



The Einstein tensor has the non-vanishing components 
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Gihih = 



I sin 2 9 f ^ r 2 m^ r + rm^ r + m fi s _ ^m^y _ | rm, r j 2 j + ^^^ 2 + ^ r ^ r ^ 2 



G — — (-2—!^- + rl ' r l,e + Tm ' r l ' e + rl ' r m ' e + 2 m ' e - 2 r ^' r + 2 ( rm,r - — ) c ; 9 
r 4r \ I 11 ml I m m / / \ m I 



With these results at hand we now derive the equations of motion. 



(91) 



B. Matter equations 

The equations of motion for the gauge field functions and for the dilaton function are obtained from the variation 
of the matter part of the action (|l|), 

= (r 2 H hr , r + H he ,e + H 2fi - n 2 j (rH 4 , r H 3 - rH 3 , r H 4 + (.fff + H\ - l) ffi)) sin 2 9 
rH 2 , r + H 1<e ~n 2 — {2HxH 3 + rH^ r )J sin 9 cos 9 - n 2 — #1 



' 9 [Ht,o + rH 2 , r ] In ( E— ] + sin 2 9 [rH hr - H 2 ,e]r\n ( E — 



m 



= (-{r 2 H 2 ^ r + H 2fifi - H he ) + n 2 j (H 4 H 3 , e - H 3 H 4fi + (H 2 + H\ - l) H 2 ) ) si,.- 
rH 1>r - ff 2 , e + n 2 — (2H 2 H 3 - H 4 . e )j sin cos + n 2 - (£T a - H 4 ) 

+ sin 2 [ril 1>r - H 2fi )\ In ( E- ] - sin 2 9 [H hB + rH 2 , r ] r In ( E— 



= (r 2 /f 3 ,r,r + ff 3 ,e,e - #3 (H 2 + H 2 ) + H x (H 4 - 2rH i>r ) - ff 4 (rfli, r - fr 2 , 9 ) + 2iT 2 if 4>e ) sin 2 
+ (F 3) e + H 2 H 4 - H\ - Hi) sin 9 cos 9- H 3 

+ sin 2 6> [# 3 , e - 1 + i? 2 i7 4 + cot 9H 3 ] In (E) fl + sin 2 9 [rH 3 , r - H^i] r In (E) 



(92) 



(93) 



(94) 



= {-r 2 H 4 ^ r - H 4 ,g : e + H 4 (H 2 + H\) + fftflg - ff 2 (1 - 2H 3 , e ) - 2rH 3 , r H 1 - H 3 {rH 1>r - H 2>e )) sin 2 9 
+ {H 2 H 3 +H X - rH hr + H 2fi - H 4fi ) sin 9 cos 9 - (H 2 - H 4 ) 

- sin 2 9 [Hi,e - H 3 H 2 - {H 2 - H 4 ) cot 9] In (E) g - sin 2 9 [rH^ r + H 4 H 3 + cot 9H 4 ] r In (E) r 

(95) 

= *,r,r + " «e 2K * f J-F 2 g + I \A F 2 r4> + \f 2 h 

r z \r z m r 2 am 91 r z v 

+ In(sin6»\/Ir 2 ) ir $ ir + ^\i\{ain9\flr 2 ) : e^,e 

(96) 

with E = e 2K *"^- The differential equations for the gauge field functions and the dilaton function are diagonal in the 
second derivatives with respect to the variable r. 

C. Metric equations 

The equations for the metric functions are the Einstein equations 

(97) 
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They represent five differential equations for three metric functions. Two of these equations are redundant and can 
be used to achieve a convenient form of the three equations to be solved numerically (see Appendix D). 

We diagonalize the Einstein equations with respect to the set of derivatives (f, r , r > m ,r,r, l,r,r, 1,6,6, l,r,e)- Except for 
the (00) equation, these equations are already diagonal. The resulting equation for the (00) component is 

1 ( r 2 f, r ,r + f,e,8 + 2rf, r ( ( ffi\ 2 ( r f,r\ \ . , a f,e lfrf >r rl, r 



(98) 

Then we solve the (it) component for (^-y 11 ) 2 and substitute this expression into the (4xj>) and (99) components, 
yielding 

8ttG 



t ( r - +T t)-^ ( 2 ( "' m " + r + '"' M -BY- (^) 2 ) + 2 ¥ + 2 ( V 2 



/ 

hi] a. 2 { — + rm ' r \ + rm > r rl ' r - Hhihi _ 2 ( Tji _ <i — I coi i 
I I \ I ml ml ml \ m I 



,m /„ 1 ( n r 2 l r , r + lga fl.e\ 2 I rl r \ 2 rl, r J, 



8. Gy (T; + T e »)^ -^f + 6^ + 4 f cot< 



(99) 



(100) 



The differential equations (|92|)-(|96|) and (|98|)-(10C) are then diagonal in the second derivatives with respect to the 
radial variable r. 

Finally we change to dimensionless quantities x, <p and 7, (eqs. (p4[)-(p6|)) and furthermore to the radial variable x, 
(eq. (^3j)), thus mapping the infinite interval of the variable x onto the finite interval [0, 1] of the variable x. For the 
derivatives this leads to the substitutions 

rF r — ► 5(1 - x)F t x , 
r 2 F^ r — » x 2 ((1 - ,ri 2 /•;,,, - 2(1 - i)F 2 ) (101) 

(102) 

for any function F in the differential equations. In this form we have solved the system of differential equations 
numerically. 

VIII. APPENDIX C 

Here we present the expansion of the functions at the origin in powers of x and at infinity in powers of 1/x. 

A. Expansion at the origin 

Expanding the gauge field functions, the dilaton function and the metric functions at the origin up to the second 
order in x we obtain 

H ± = x 2 sin 6 cos 9H n + 0(x 3 ) , (103) 

H 2 = l-x 2 ^(l~2sin 2 9) + 0(x 3 ) , 
H 3 = x 2 sin 9 cos 9H 31 + 0(x 3 ) , 

H 4 = 1 - x 2 - H 31 sin 2 ej + 0(x 3 ) , 

V = ^-x 2 f(l-3cos 2 9)^~n 2 1S in 2 9^-e 2 ^(H 11 -2H 31 ) 2 ^j + 0(x 3 ) , 
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/ = ./'„ ( I - x 2 ( A (i _ 3 cos 2 0) - n 2 sin 2 6-^-e 2 ^° (H n - 2H 31 ) 2 ) j + ()(./' i . 



/// = /„ ( I x 2 ih (l-2cos 2 6») - n 2 sin 2 6»^e 2wo (H n - 2H 31 ) 2 ) ) + ()(.r ; ) 



l = l I 1 -® J (l-4cos^6»)J + 0(a; d ) , 
where i?n, #31, tfo, ifi, fa, /1, and l\ are constants. 

B. Expansion at infinity 

Expanding the gauge field functions, the dilaton function and the metric functions at infinity in powers of ^ we 
obtain 

Hi = -^H 12 sin 9 cos 6 + O (J^ J , (104) 
H2=±1 + (cos 2 9 - sin 2 9) + O (i) , 

i? 3 = isin6»cos6l J ff3i + ^sin6lcos6l(±2iJi 2 + iJ3i (/1 + 27^1)) + O 

Hi=±(l+ ifl- 31 sin 2 - JL sin 2 (±2i? 12 + H 31 (/1 + 2 7 ^))^ + O ^ 
1 „ / 1 



y = — >fii + O 



1 - 1 9 / 1 

m = 1 + -^-Z 2 + -r sin 2 9m 2 + O 

x z x 

1 = 1 + \i 2 + o ' J 

with constants -ffi2, Sai, /1, <£>i, ^2 and l 2 . 

IX. APPENDIX D 

Here we discuss some technical aspects concerning the construction of the axially symmetric solutions. 

A. Metric 

Prior to the metric (^) we have employed the isotropic metric in the form 1 30 1 

ds 2 = -fdt 2 + (dr 2 + r 2 d9 2 ) + lr 2 sin 2 9d<p 2 , (105) 
where the metric functions /, /I and I are functions of the coordinates r and 9. At infinity, the functions satisfy 

/ = 1 + -fi + O ( \ 

x \x 

fi = --h + O f ' 

x 
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l = l_I/ 1+0 (_L) , CIDli) 



x \ x , 

(107) 

with constant f\. We have also replaced the function p, by the new function to, 

m = e' 2 . (108) 

Our final choice of metric parametrization with 

to = y and f = j , (109) 

has the advantage that it does lead to a faster convergence of the metric functions to and I to their limiting values 
(see Appendix C). 

B. Einstein equations 

We have diagonalized the Einstein equations with respect to the second derivatives in r of the metric functions. But 
this diagonalization is ambiguous, since only 3 of 5 equations are needed. In the earliest calculations we followed the 
diagonalization scheme employed in ref. However, these calculations did not properly produce the asymptotic 

behaviour of the derivative d x ji. Using another linear combination of the equations, the problem was shifted to the 
asymptotic behaviour of the derivative d x f. The problem persisted even for the spherical solutions, indicating that 
its origin resided in the program package FIDISOL. Therefore we looked for a combination of the equations which 
FIDISOL would deal with properly. We realized, that elimin ating the square of the first derivatives of /, {d r f) 2 , from 
the metric equations, would precisely do this. Eqs. ( p9j ) and ( |100[ ) then become identical in the spherically symmetric 
case, n = l. In this form the computations with FIDISOL give the correct behaviour of the functions for n = 1 and 



also for 7i > 1. This motivates our final choice of the equations, (B8f)-100) 



C. Regularity condition 

Regularity on the z-axis requires the condition (||) to be satisfied. When the numerical calculations do not impose 
this condition, it is not exactly satisfied. Typical values for the ratio m/l are 0.995 instead of precisely one. To strictly 
satisfy the regularity condition (||) we have introduced a new function 

TO 

9= j, (HO) 

imposing condition (^) as boundary condition for this function. The resulting solutions respect the regularity condition 
and are otherwise in excellent agreement with the solutions obtained before. 



22 





A* 


-j-J— 

k / n 


1 


2 


3 


4 


i 
i 


0.577 


0.961 


1.297 


1.607 


2 


0.685 


1.262 


1.770 


2.239 


Q 
O 


0.703 


1.365 


1.976 


2.549 


4 


0.707 


1.399 


2.063 


2.698 




0.707 


1.414 


2.121 


2.828 




tmax (pmax) 


k/n 


l 


2 


3 


4 


1 


1.075 (0.) 


0.177 (0.90) 


0.098 (1.59) 


0.072 (2.37) 


2 


11.63 (0.) 


0.910 (0.30) 


0.380 (0.66) 


0.235 (1.10) 


3 


79.70 (0.) 


3.443 (0.09) 


1.124 (0.28) 


0.601 (0.51) 


4 


498.2 (0.) 


12.01 (0.03) 


3.064 (0.11) 


1.435 (0.24) 



Table 1 

The dimensionless mass /Lt, the maximum of the energy density e max and the location p max of the maximum of the 
energy density of the EYMD solutions of the sequences n = 1 — 4 with node numbers k = 1 — 4 and 7=1. For 
comparison, the mass of the limiting solutions is also shown. Note, that fi/n decreases with n for fixed finite k. 
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A* 


-j-J— 

k / n 


1 


2 


3 


4 


i 
i 


.829 


1.385 


1.870 


2.319 


2 


.971 


1.796 


2.527 


3.197 


Q 
O 


.995 


1.935 


2.805 


3.621 


4 


.999 


1.980 


2.921 


3.824 




1. 


2. 


3. 


4. 




tmax (Pmax) 


k/n 


l 


2 


3 


4 


1 


1.23 (0.) 


0.26 (0.7) 


0.16 (1.3) 


0.12 (1.8) 


2 


2.55 (0.) 


0.49 (0.2) 


0.26 (0.5) 


0.19 (0.8) 


3 


2.92 (0.) 


0.62 (0.07) 


0.35 (0.2) 


0.25 (0.4) 


4 


2.98 (0.) 


0.68 (0.02) 


0.40 (0.1) 


0.28 (0.2) 



Table 2 

The dimcnsionless mass /i, the maximum of the energy density e max and the location p max of the maximum of the 
energy density of the EYM solutions of the sequences n = 1 — 4 with node numbers k = 1 — 4. For comparison, the 
mass of the limiting solutions is also shown. Note, that fj,/n decreases with n for fixed finite k. 
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EYM 


EYMD 


YMD 


k/f 





0.5 


1.0 


2.0 


oo 


1 


1.870 


1.659 


1.297 


0.811 


1.800 


2 


2.524 


2.250 


1.770 


1.114 


2.482 


3 


2.805 


2.505 


1.976 


1.247 


2.785 


4 


2.922 


2.611 


2.063 


1.304 


2.913 




3 


2.683 


2.121 


1.342 


3 



Table 3 

The dimensionless mass // of the EYMD solutions with winding number n — 3 and up to 4 nodes for several values of 
the dilaton coupling constant 7. For comparison, the last row gives the mass of the limiting solutions, the first column 
gives the mass of the EYM solutions, and the last column the scaled mass of the corresponding YMD solutions 
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n = 1 7=1 




-4-2 2 

log 10 (x) 



Fig. 1 

The coordinate transformation between the isotropic coordinate x and the Schwarzschild-like coordinate x is shown 
for the spherically symmetric solutions (n = 1) of EYMD theory with 7 = 1 and k = 1 — 4. Also shown is the 

coordinate transformation for the limiting EMD solution. 
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n = 1 7=1 




Fig. 2a 

The metric function / is shown for the spherically symmetric solutions (n = 1) of EYMD theory with 7=1 and 
k = 1 — 4. Also shown is the metric function of the limiting EMD solution. 
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Fig. 2b 

Same as Fig. 2a for the metric function m. 
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Fig. 3a 

Same as Fig. 2a for the gauge field function w. 
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n = 1 7=1 
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Fig. 3b 

Same as Fig. 2a for the dilaton function ip. 
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n = 1 7 = 




-4-2 2 

log 10 (x) 



Fig. 4 

The coordinate transformation between the isotropic coordinate x and the Schwarzschild-like coordinate x is shown 
for the spherically symmetric solutions (n = 1) of EYM theory with k = 1 — 4. Also shown is the coordinate 

transformation for the limiting RN solution. 
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n= 1 7=0 




-4-2 2 

log 10 (x) 



Fig. 5a 

The metric function / is shown for the spherically symmetric solutions (n = 1) of EYM theory with k = 1 — 4. Also 

shown is the metric function of the limiting RN solution. 



32 



n= 1 7=0 



/ / 



/ / / 
/ / / 

/ 



/ /• 



/ / / k = 2 



/ / 1 

i 1 : k=3 

. — 



k = 4 



./ 

/ / 



-4-2 2 

log 10 (x) 

Fig. 5b 

Same as Fig. 5a for the metric function m. 
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f(p ; z) n=2 k=l 7=1 




Fig. 6a 

The metric function / of the EYMD solution with 7=l,n = 2,fc = lis shown as a function of the compactified 

coordinates p and z. 
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l(p,z) n=2 k=l 7=1 




p 



Fig. 6c 

Same as Fig. 6a for the metric function 
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H s (p,z) n= S k=l 7 = 1 
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H g (p,z) n= S k=l 7 = 1 




39 



H 4 (p,z) n= S k=l 7 = 1 
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Nodes n = 2 k=l 7= 1 
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42 



e(p,z) n=2 k=l y = l 
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e(p,z) n=2 k=l 7=0 




Fig. 10 

The energy density of the matter fields e of the EYM solution with n = 2, k = 1 is shown as a function of the 

compactificd coordinates p and z. 
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Fig. 11a 

The metric function / of the EYMD solutions with 7 = 1, n = 1,2, and 4 and fc = 1 is shown as a function of the 

coordinate x for the angles 9 = 0, 7r/4 and ir/2. 
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Fig. lib 

Same as Fig. 11a for the metric function to. 
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Fig. 11c 

Same as Fig. 11a for the metric function I. 
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Fig. 12b 

Same as Fig. 11a for the gauge field function H 2 - 
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Fig. 12c 

Same as Fig. 11a for the gauge field function H 3 . 
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Fig. 12d 

Same as Fig. 11a for the gauge field function H 4 . 
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Fig. 13 

Same as Fig. 11a for the dilaton function 
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Fig. 14 

Same as Fig. 11a for the energy density of the matter fields e. 
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Fig. 15a 

The metric function / of the EYM solutions with n = 1,2, and 4 and k = 1 is shown as a function of the coordinate 

x for the angles 9 = 0, 7r/4 and ir/2. 
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Fig. 15b 

Same as Fig. 15a for the metric function to. 
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Fig. 15c 

Same as Fig. 15a for the metric function I. 
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Same as Fig. 15a for the gauge field function H 2 - 
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Same as Fig. 15a for the gauge field function H 3 . 
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Fig. 16d 

Same as Fig. 15a for the gauge field function H 4 . 
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Same as Fig. 15a for the energy density of the matter fields e. 
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Fig. 18a 

The metric function / of the EYMD solutions with 7 = 1, n = 2, and fc = 1 — 4 is shown as a function of the 
coordinate x for the angles 9 = 0, 7r/4 and 7r/2. Also shown is the metric function of the limiting EMD solution. 
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Same as Fig. 18a for the metric function to. 
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Fig. 18c 

Same as Fig. 18a for the metric function /. 
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Fig. 19a 

Same as Fig. 18a for the gauge field function Hi. 
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Same as Fig. 18a for the gauge field function H 2 - 
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Fig. 19c 

Same as Fig. 18a for the gauge field function H 3 . 
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Same as Fig. 18a for the gauge field function H4. 
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Same as Fig. 18a for the dilaton function ip. 
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Fig. 22a 

The metric function / of the EYM solutions with n = 2, and k = 1 — 4 is shown as a function of the coordinate x for 
the angles 9 = 0, 7r/4 and 7r/2. Also shown is the metric function of the limiting RN solution. 
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Fig. 22b 

Same as Fig. 22a for the metric function m. 
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Fig. 22c 

Same as Fig. 22a for the metric function I. 
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Fig. 23 

Same as Fig. 22a for the energy density of the matter fields e. 
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Fig. 24a 

The logarithm of the absolute deviation from the limiting solution Afik = Moo(oo) — /Zfc(oo) for the EYMD masses as 
a function of the node number k with dilaton coupling constant 7 = 1 for winding number n = 1 — 4. The straight 

lines are fitted to the data with k = 3 and k = 4. 
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Fig. 24b 

Same as Fig. 24a for the EYM solutions, i. c. 7 = 0. 
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